Astronomy &: Astrophysics manuscript no. nis0137 


February 2, 2008 


(DOI: will be inserted by hand later) 





On the thermal behaviour of small iron grains 
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Abstract. The optical properties of small spherical iron grains are derived using a Kramers-Kronig-consistent 
model of the dielectric function including its dependence on temperature and size. Especially discussed is the effect 
of the size dependence, which results from the limitation of the free path of the free electrons in the metal by 
the size of the grain, on the absorption behaviour of small iron spheres and spheroids. The estimated absorption 
properties are applied to study the temperature behaviour of spherical and spheroidal grains which are heated by 
the interstellar radiation field. 
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1. Introduction 

According to the cosmic abundances given by Anders 
& Grevesse ipg^ . iron is the ninth most abundant 
element in the interstellar medium (ISM) and should 
therefore be an essential constituent of interstellar dust 
grains. Iron particles are expected to form independently 
of the oxygen-to-carbon ratio in circumstellar shells (see 
Jones 119901 and references therein) and have been sup- 
posed to be, apart from silicates, the main condensates in 
the metal-rich core region of supernova ejecta (Hoyle & 
Wickramasinghe Fl 970|l . SN 1987A was the first supernova 
where the condensation of new grains in ejected material 
could be observed. The composition of these grains could 
not be clarified but indications, summarised by Wooden 
H1997|l . point to predominantly iron-rich grains. 

Even though the ISM will probably be enriched in 
iron particles it is uncertain whether they can survive 
and play an important role as a dust component of in- 
terstellar grains (Dulcy 1980, Jones 1990). If they are 
present, iron particles may affect the interstellar extinc- 
tion (Wickramasinghe & Nandy 11971(1 , the chemistry in 
the ISM due to catalytic processes on grain surfaces 
(Tabak I1978|l . and may contribute to the diffuse emis- 
sion from interstellar dust grains (Chlewicki & Laureijs 
I1988|l . Further, iron grains may be an important grain 
species to explain the Rosseland mean opacity in molec- 
ular clouds cores (above « 575 K) and accretion disks 
(above « 425 K) (Pollack et al. lTM^i . On the other hand 
it has been claimed, based on observations at 90 GHz, that 
not more than 5% of interstellar iron is in form of metallic 
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iron grains or inclusions (Draine & Lazarian. ll999|l . This 
would support the hypothesis that the survival time of 
a relatively pure iron grain is quite short. However, even 
though overall iron may be not an important constituent 
of interstellar grains it does not exclude the possibility of 
local enhancements where these grains may form. 

Mathis, Rumpl & Nordsieck 11977|l found, from analy- 
sis of the extinction curve, that the interstellar medium is 
mainly composed of grains much smaller than ~ 0.25 /xm. 
Very small grains are furthermore thought to be respon- 
sible for the diffuse IR emission at wavelengths shorter 
than ^ 60 /im (see e.g. Draine & Anderson 11985(1 due 
to the stochastic heating process, where very small grains 
reach temperatures much higher than the corresponding 
equilibrium temperature. To model their emission or their 
ir optical properties, it is usually assumed that the dielec- 
tric function derived for solids can also be applied for very 
small grains. At least for small metallic grains, as shown 
by experiments with very small silver spheres, it is known 
that the dielectric function has to be corrected for grain 
size (Kreibig & Fragstein 1969. Kreibig J974J . The influ- 
ence of the grain size on the dielectric function can be even 
stronger than its dependence on temperature. One expects 
that the optical properties of small metallic spheres should 
be much less affected by temperature than those of larger 
grains as shown, for example, by calculations of the ab- 
sorption behaviour of spherical graphite grains (Draine & 
Lee lTOMj) . 

Here, we present calculations for the temperature be- 
haviour of pure iron grains in the ISM heated by the in- 
terstellar radiation field (ISRF) where we have taken into 
account both the temperature and the size dependence 
of the dielectric function. It has been suggested that a 
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part of the interstellar grains may have a non-spherical 
shape (see Voshchinnikov et al. 1999). Therefore, apart 
from simple iron spheres, spheroidal iron grains will also 
be considered. The influence of the shape of small grains 
of different compositions, including iron, on temperature 
has also been discussed by Voshchinnokov et al. (^999). 
Their results are based on the dielectric function of solids 
using a sophisticated code to derive the optical properties 
of spheroidal grains. Here a much simpler approach is used 
demonstrating the effect of the size dependence of the di- 
electric function on the cooling behaviour of spheroidal 
iron grains. 

In Sect. H we describe the model used to derive the 
optical properties of small iron spheres and analyse in 
more detail the absorption behaviour of small spherical 
and spheroidal iron grains. In Sect. 13 the optical proper- 
ties are used to derive the temperature behaviour of spher- 
ical and spheroidal small iron grains in the ISM which are 
heated by the ISRF. The derived temperature distribu- 
tions are used to obtain the emission spectrum of spherical 
iron grains assuming a grain size distribution as proposed 
for the ISM including the stochastic heating process of 
very small grains. The results are discussed in Sect. 4 and 
summarised in Sect. 5. 

2. The optical properties of iron grains 

The interaction of an electromagnetic wave of a given en- 
ergy with a grain is generally determined by the response 
of the corresponding solid, expressed by the complex di- 
electric function e, the shape, and the size of the grain. 

As for non- metallic grains the magnetic permeability is 
taken to be /i = 1. As will also be shown this assumption 
should be sufficient to describe the cooling of iron grains in 
the ISM, which is the main purpose of this paper. The con- 
tribution of the magnetizability of metallic grains to the 
optical properties has been studied by Draine & Lazarian 
l|1999|l and is only significant at wavelengths larger than 
~ 1000 fi-ai. 

2.1. A model for the dielectric function 

To derive the optical properties of spherical iron grains 
we used a self consistent simple model for the complex 
dielectric function e{uj) = ii{io) + ie2(^) based on optical 
data derived for bulk iron at a temperature of 298 K. 

We followed basically the method used by Draine & 
Lee {!1984|) to determine the optical properties of small 
graphite grains and separated the dielectric function in 
Se-f and Se'', which correspond to the free and the bound 
electrons of the grain: 

e^l + Sef+Se\ (1) 

It is assumed that the dependence on temperature and 
grain size is given through the free electron part which is 
described by a single Drude term 



where F is the collision rate and cOp the plasma frequency. 

The parameters of the Drude term of bulk iron are 
determined at long wavelengths where the free electron 
gas dominates the optical behaviour of metals. To derive 
a self consistent dielectric function we have chosen the 
imaginary part as given by £3 = £2 — and used the 
Kramers-Kronig relation (Landau & Lifschitz, I1985|l to 
estimate the real part e\ . 

The main dependence of 6e-^ on size and temperature 
is given by the collision rate F which can be taken as the 
sum of the collision rates Fj and Fg due to scattering events 
in the solid and at the outer surface of the grain (Kreibig 
& Fragstein 119691 Kreibig fW^ . A further temperature 
dependence is assumed to be caused by the thermal ex- 
pansion of the solid which alters the electron density Uc 
and therefore LUp (wp cx nl^'^). 

Inside the solid, electrons are scattered at impuri- 
ties, defects and phonons, but only the phonon scattering 
events lead to a strong temperature dependence. In the 
case of bulk iron, the scattering at impurities and defects 
is negligible with the exception of very low temperatures 
where they are the dominant scattering process inside the 
grain. 

As a first approximation the scattering rate at phonons 
may be assumed to be independent of the amount of impu- 
rities and defects and the scattering rate at impurities and 
defects independent on temperature so that the effects can 
be easily de-coupled. Under these assumptions the colli- 
sion rate due to impurities and defects can be scaled by 
a factor ( to allow consideration of the cases where the 
iron grains are either less pure or have more defects in 
comparison to pure iron grains. 
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Fig. 1. Critical grain sizes of iron spheres where the col- 
lision rate Fi of the free electrons in the material is equal 
to the collision rate at the surface shown for C = 1, 10, 
and 100. Below the curves the collision rate is dominated 
by surface scattering. 

The total collision rate is given by 
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where (Jo{T, (^) is the conductivity of the iron grains taken 
to be 



ao(r, C) a\ ab"ik(T) a^-^^iO K) J ' ' ' 

(Tg""^(r) is the adopted temperature-dependent conduc- 
tivity of bulk iron and a = CTo(298 K)/ct^"1''(298 K) a 
(smaU) correction factor where cro(298 K) is the conduc- 
tivity obtained by fitting the Drude term to the optical 
properties of bulk iron (see below), t;^ is the Fermi ve- 
locity of the electrons and a the radius of the grain. The 
parameter f3 is of order one and is given by the scattering 
process at the surface. As Draine & Lee we have taken 
P = 1, which is equivalent to isotropic scattering. 

The collision rate at the surface becomes important 
when the free path of the electrons is greater than the size 
of the grain. Where the collision rate in the solid is equal 
to the collision rate at the surface is shown in Fig. for 
different values ^. In the case of pure very cold iron grains 

= 1) the scattering should be dominated by surface 
scattering for all grains with radii smaller than ^ 20 /im. 
If the grains are less pure and have more defects the effect 
is limited to smaller grain sizes. 

As the purity of 'astrophysical' iron grains is not 
known we considered pure iron grains. The collision rate 
assumed in this paper can therefore be regarded as a lower 
limit of the collision rate of actual iron grains. However, 
this assumption should not affect interstellar iron grains 
which, if they do exist, are believed to be quite small with 
maximum sizes not much larger than ~ 0.1 /im. As seen in 
Fig. n] their mean free path should be determined by the 
size of the grain unless the grains are highly disordered or 
impure with ( > 100. The effect of the purity and defects 
on the temperature behaviour of larger iron grains will be 
discussed in Sect. 14.1.11 

The optical data used to derive the optical properties 
can be taken from Table ^ We interpolated data from 
Lynch & Hunter H199f)|l and Moravec et al. H197()|l in the 
energy range 51.57 eV to 26 eV using f{E) = cE^ . A 
smooth transition was made in the range 3.5 to 5 eV be- 
tween data from Moravec et al. and the data published in 
the 'CRC Handbook for Chemistry and Physics' (1191)81). 
Above 0.8 eV we used the averaged data from Ordal et al. 



Table 1. Used optical data 



energy intervall [eV] 


reference 


10 000 


- 51.57 


Lynch & Hunter (.1996,) 


51.57 


- 26.0 


interpolation (see text) 


26. 


- 5 


Moravec et al. I|1976|l 


3.5 - 


' 1.55 


CRC (I998| 


1.55 - 


0.0124 


Ordal et al. (IM^ 


0.0124 - 


' 0.00124 


extrapolation (see text) 



From the optical data in the range of 1 to 100 fim we 
determined the plasma frequency tOp and the internal col- 
lision rate Fj of bulk iron at room temperature which are 



given in Table [3 To prevent amplification of the electro- 
magnetic wave inside the solid we have taken care that 
£2 = £2 — ^2 > in the considered energy range. The de- 
rived parameters of the Drude term (up = 27 020 cm^^, 
Fj — 125.2 cni^^) are slightly smaller than the ones 
used by Ordal et al. for extrapolation above 100 fim 
{ijjp = 29 500 cm~^ and Fi = 156 cm~^) but very close 
to those used by Draine & Lazarian H1999|l . 

The DC conductivity of iron at room temperature is 
found to be ctq = cjp/47rFi = 7.65 • 10^'' s~^ which is only 
about a factor of 1.2 lower than the measured value of 
9.12 • lO^s s-i (CRCHnnni)- For comparison in the case of 
graphite the DC conductivity given by the parameters of 
the Drude term is a factor of 3 smaller than the measured 
value (Draine & Lee 11984(1 . 



Table 2. Parameters used for the dielectric function (5e' 



parameter 


value 


plasma frequency (298K)'' 


iOp [s ^] 


5.090 • 10''' 


internal collision rate (298K)" 




2.693 • 10^^ 


Fermi velocity'' 


VF [cm/s] 


1.98 ■ 10* 



" Fit parameters with single Drude term. 
' Taken from Ashcroft & Mermin llT^T^ . 



The derived dielectric function is shown in Fig. [51 and 
Fig. 131 At long wavelengths (Fig. [SJ the functional form 
of is qualitatively similar to the dielectric function e\ 
of graphite perpendicular to the main axis (Draine & Lee 
^84), which shows metal-like behaviour. The imaginary 
part e\ increases with wavelength and reaches its highest 
values at ~ 60 /im. Beyond 100 /xm we assumed that the 
values of e\ decrease with a power law. The explanation of 
the behaviour of the dielectric function in the IR is out 
of the scope of this paper and it cannot be excluded that 
this behaviour is a result of the simplification in modelling 
the free electron regime with a single Drude term. 

The temperature dependent DC conductivities 
cr;j'""'(r) of bulk iron are taken from the 'CRC Handbook 
of Chemistry and Physics' H1998|l . The expansion coeffi- 
cients of iron for correcting the plasma frequency with 
temperature was taken from the 'American Institute of 
Physics Handbook' H1972|l . For the collision rate at the 
surface we used for the free electron gas of iron a Fermi 
velocity vp = 1.98 • 10^ cm/s (Ashcroft & Mermin [WTSjl . 

2.2. Optical properties of spherical pure iron grains 

We estimated the optical properties of spheres above 
27ra/A = 0.1 using the MIE-formalism from Bohren & 
Huffman 1(1 983(1 . At longer wavelengths we used the dipole 
approximation for spheres presented by Draine & Lee 
(^.984). The results for iron grains at a temperature of 
25 K are shown in Fig. 01 At this temperature the mean 
free path of the electrons would be larger than the grains 
so that the collision rate F is given by the collision rate at 
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Fig. 2. Real and imaginary part of the dielectric function for individual pure small iron spheres at energies above 
1 eV. For comparison also the real and imaginary parts of e*^ are shown as thin solid lines. The part of 63 shown as thin 
dotted line are extrapolated values used for the Kramers-Kronig relation to derive e^". The real part of the dielectric 
function is always negative except for the three branches shown. The short and the long dashed thin lines in the right 
hand figure refer to the dielectric function 6el — Im(e — e'') of an iron sphere with a = 0.001 /im and a — 1.0 ^m, 
respectively. 
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Fig. 3. Real and imaginary part, given as conductivity a — e2w/47r, of the dielectric function e above leV at a 
temperature of 298 K. The chosen values of the imaginary part and the result for the real part after Kramers-Kronig 
analysis are shown as solid lines. The values differ somewhat from the published values. The derived values of the real 
part of the dielectric function lie between those from Johnson & Cristy 1)1974(1 and those from Moravec et al. H1976|l . 
The real part e\ = Re(e — 5e^) (shown as dashed-dotted line) is the direct result derived with the Kramers-Kronig 
relation from the imaginary part £2 — ^2 ~ ^^2- 



the grain surface and therefore through T w vp/a. While 
scattering (Qscaj 5sca) is almost insensitive to the grain 
temperature, the absorption of larger iron grains can have 
large temperature variations. 

The dependence of the absorption behaviour of iron 
grains on temperature is shown in Fig. [S] for two different 
grain sizes with radii of 0.01 and 1 /im. For comparison 
we also show the absorption cross section (C^^) due to 
magnetisation effects inside iron spheres well below the 
Curie temperature as derived by Draine & LazarianUHMI 
If correct this effect clearly dominates the absorption of 



cold iron grains at wavelengths longer than ~ 1000 /im. 
Above the Curie temperature the iron spheres become non 
magnetic and the absorption at long wavelength should be 
determined by the electronic response of the material. As 
the magnetisation does not affect the cooling behaviour of 
the grains this effect will be ignored in the following. 

Due to the high Fermi velocity only the optical proper- 
ties of larger iron grains with a ^ 0.01/im are affected by 
grain temperature. Even a variation of grain temperature 
over a huge range (10 - 1200 K) causes only a slight change 
in the absorption behaviour of grains with a = 0.01 /xm. 
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Fig. 4. Derived optical properties of spherical iron grains. 
Shown are the absorption coefRcients (5abs(A,a), scatter- 
ing coefficients Qsca{X, a), and the scattering parameter 
(7sca(A, a) = {cos{d)) at 25 K for various grain radii from 
10 A to 10 /im with a difference in radius of A log a[/^m] — 
0.5. Thick lines are labeled with the corresponding grain 
radii in microns. 
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Fig. 5. Effect of the temperature on the absorption be- 
haviour of spherical iron grains. For various temperatures, 
the absorption coefficients of iron grains with radii 1.0 and 
0.01 /Ltm for n — 1 (thick lines) are shown. For comparison 
the absorption coefficients due to magnetisation effects for 
magnetic single and multi-domain iron grains well below 
the Curie temperature (744° K) are given (thin dashed 
and solid line). 



Grains with a radius of a = 1 /im on the other hand have, 
at a wavelength of 100 /im, an order of 10 higher absorp- 
tion at 300 K than at 10 K. In this wavelength region 
the absorption coefRcients increase with temperature. At 
longer wavelengths it is possible that colder iron grains 
emit more readily than warm grains do. 



2.2.1. Dipole absorption of spherical pure iron grains 

As is well known, metallic spheres show relatively com- 
plex absorption characteristics at long wavelengths. The 
behaviour is in contrast to non-metal grains where the 
IR-absorption is close to Qabs oc a/)?. This is due to the 
fact that at long wavelengths the absorption behaviour of 
grains with a high dielectric function, as for metals, can 
be dominated by the magnetic dipole absorption caused 
by eddy currents which are created through the changing 
magnetic field inside the conductor (Landau & Lifschitz 
11985(1 . In our calculation the absorption behaviour of iron 
spheres at wavelengths shorter than 60 /zm is addition- 
ally partly affected by the high values of the dielectric 
function e''. 

However, as seen in Fig. 0] and Fig. |S1 inside certain 
regions of the parameter space defined by the wavelength 
A and the grain size a the absorption coefficient can still 
be described by a simple approximation Qabs(a, A) (x a^'A* 
with fixed values for s and t. The absorption properties at 
long wavelengths can be explained using the electric and 
the magnetic dipole approximation for the cross sections 
Cabs = T^o?Qahs and are discussed in some detail in the 
Appendices El and El 

2.3. Dipole absorption of small spheroidal pure iron 
grains 

As a special case of non-spherical iron particles we will 
consider spheroidal iron grains. In adopting the model 
of the dielectric function described in Sect. 12.11 to non- 
spherical grains one has to choose the appropriate value 
of the collision rate in Eq. |21 As an approximation of the 
mean free path due to surface scattering we have taken 
the averaged distance from the center to the surface: 

, , b ( arcsin(e) a > b (prolate), , . 

^""^ " ; \ In (e -f V^^TT) b > a (oblate), ^ ' 

where = |1 — (&/a)^|. a and b are the semi-major axis 
parallel and vertical to the rotation axis. For a sphere 
(r) is equal to the mean free path of isotropic scattering 
(/3 = 1). In the limit of an elongated prolate grain (e « 1) 
the averaged distance from the centre to the surface is 
equal to 7r6/2 and only slightly different to the mean free 
path of 26 in the case of an infinite needle if isotropic 
scattering is assumed. 

In general for the calculation of the optical proper- 
ties of non-spherical grains, even for spheroidal grains, can 
be very complicated (see e.g. Voshchinnikov & Farafonov, 
11993(1 . Here we are interested in the temperature be- 
haviour of grains in the interstellar medium where the ra- 
diation is almost isotropic and unpolarised. In such cases 
it is sufficient to use only the cross section averaged over 
the orientations of the spheroids. 

For simplicity, we considered spheroidal grains small in 
comparison to the wavelength of the absorbed photon so 
that the optical properties can be described by the dipole 
approximation. For the electric dipole absorption we used 
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Fig. 6. Angle-averaged absorption cross section (Cabs) of prolate (a > b) and oblate (a < b) iron grains with volumes 
equivalent to the volume of a sphere with radius — 0.01 /im. For comparison also the absorption of a sphere is given. 
Solid and dashed lines correspond to absorption where the size and temperature dependence of the dielectric function 
has been taken into account and has been ignored, respectively. The dotted line shows the electric dipole absorption 
of a sphere for the case of temperature and size independent dielectric function. 



the approximation for ellipsoidal grains and for the mag- 
netic dipole absorption we adopted the expression given 
for spheres (Draine & Lee. 1198 jjl . 

The angle averaged absorption cross sections of 
spheroidal iron grains for a number of different ratios 
a/b are shown in Fig. |H| (solid lines). The volume of the 
spheroids is chosen to be that of a sphere with a radius 
of r„ = 0.01 iim. For comparison also the cross sections 
for temperature and size independent dielectric function 
(dashed lines) is given. 

The absorption behaviour at long wavelengths of a 
spherical iron grain with a = 0.01 ^m is dominated by the 
electric dipole absorption if the size dependence is taken 
into account. This is also true for the spheroids as the con- 
tribution of the electric dipole absorption only increases 
with axis ratio while that of the magnetic dipole absorp- 
tion is assumed to be constant. 

As can be seen in the figure, by considering the size 
dependence of the dielectric function, the absorption at 
long wavelengths is clearly enhanced in comparison to 
the calculation where this dependence is ignored. On the 
other hand a resonance absorption in the IR due to sur- 
face modes (best seen for a prolate grain with a/b — 100) 
is reduced. 

The absorption properties at long wavelengths is qual- 
itatively different if the magnetic dipole absorption be- 
comes important, which is the case of larger iron grains. 
The behaviour should be similar to the case shown in 
Fig. where the size and temperature dependence of the 
dielectric function is ignored. As can be seen in the fig- 
ure, at long wavelengths the absorption behaviour of the 
sphere is then dominated by the magnetic dipole absorp- 
tion. Because of the assumption of magnetic dipole ab- 
sorption, the absorption efficiency at those wavelengths is 
almost insensitive to the axis ratio a/b unless the electric 
dipole absorption becomes important. 



The absorption behaviour in the optical on the other 
hand increases with increasing axis ratio as long as the 
ratio is small. Therefore, slightly prolate or oblate small 
iron grains with a volume corresponding to a sphere with 
ry > 0.015 /im can absorb optical light more efficiently 
than spheres but show a similar cooling behaviour. If the 
grains are heated by the ISRF for example this results in 
an increase of the grain temperature with increasing axis 
ratio fSect l3.1.2|l . 

3. Iron grains in the ISIVI 

As an example we considered spherical and small 
spheroidal iron grains which are heated by the ISRF. For 
the radiation field we have taken the energy density in 
the solar neighbourhood as derived by Mathis, Mezger & 
Panagia fMIylP lT^g^ . scaled by a factor x- 

ux^X (^^r'' + E W,^Bx{T,)^ + ^Sa(2.7 K) (6) 

with dilution factors W2 = lO^^'', W2 = 10'^^, W3 = 4 x 
10^^^ and the corresponding temperatures T2 = 7500 K, 
T3 = 4000 K and T4 = 3000 K. u^^® is the energy den- 
sity of the ultraviolet photons in the solar neighbourhood 
(MMPllMS Tab. CI). The maximum energy of the pho- 
tons is 13.6 eV. 

3.1. Grain temperature 

First we ignore the temperature variation in the case of 
very small grains and characterise all iron grains with sin- 
gle temperatures T^q where the heating is equal to the 
cooling of the grain: 

j dXcUX (Cabs) = j dA4^BA(req) (Cem) , (7) 
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where c and B\{Tcq) are the velocity of Hght and the 
Planck function. (Cabs) and (Ccm) are the absorption cross 
sections averaged over the orientations of the grain with 
(^6111) = (Cabs)- In Sect. 13.1!^ we consider the temperature 
variation of small iron grains due to stochastic heating. 

3.1.1. Equilibrium temperatures of iron spheres 

For the photon density ux we assumed three different 
strengths x = 0.1 x = ^ and x = 10- For comparison 
we also derived the equilibrium temperatures for spheres 
where we assumed a dielectric function which is indepen- 
dent of temperature and size. The derived equilibrium 
temperatures of iron grains with radii ranging from 10 A 
to 10 fim are shown in Fig.d 

As mentioned by Chlewicki & Laureijs (|1988'), the 
equilibrium temperature of iron grains in the ISM is a 
strong function of grain size. Where the magnetic dipole 
effect becomes important at a « 0.015 /im the tempera- 
ture falls from relatively high to low values which are close 
to temperatures of silicate or graphite grains. The strong 
temperature variation with grain size is also expected for 
collisionally heated iron spheres (see Fischera et al. l2002|) . 
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Fig. 7. Equilibrium temperatures of spherical iron grains 
in the ISM for three different radiation strengths (x — 0.1, 
X = 1 and X — 10). Dashed and dotted lines correspond 
to calculations in which we included and ignored the tem- 
perature and size dependence of the dielectric function. 

The grain temperatures are so low that the collision 
rate of the free electrons is dominated by surface scatter- 
ing (see Fig. . As a consequence of the absorption be- 
haviour in the IR fSect. I2'!TT|| the cooling of the smallest 
grains increases while the cooling of larger grains becomes 
less efficient. Therefore, if the size dependence is taken 
into account, the temperature of smallest iron spheres is 
lower and the temperature of larger iron spheres higher 
compared with the case where the size dependence of the 
dielectric function has been ignored. We mention that the 
cooling of the smallest grains could be even stronger if 



the dielectric function in the IR were lower than assumed 
here. 

3.1.2. Equilibrium temperatures of spheroids 

For comparison we considered spheroids with a radius 
Tv = 0.01 /xm and Tv = 0.05 /xm. The derived tempera- 
tures of these iron grains are shown in Fig. |S1 We also give 
the equilibrium temperatures of iron spheroids for a dielec- 
tric function ignoring the temperature or size dependence. 
The temperatures for extremely elongated or flattened 
iron grains, in particular for the grain with ry = 0.05 /im, 
should be taken as approximate values because of the fact 
that in this case the electric dipole approximation in the 
UV gives only a poor representation of the actual absorp- 
tion behaviour. 

Taking the size and temperature dependence of the 
dielectric function into account, the two sizes considered 
show qualitatively a very different temperature behaviour. 
We found that the temperature of a spheroid with a ra- 
dius Tv = 0.01 /xm becomes smaller with an increasing 
axis ratio. This behaviour is therefore the same found for 
spheroids of other grain compositions (Voshchinnikov et 
al. I1999|l . On the other hand the temperature of an iron 
grain with = 0.05 increases first before the temperature 
decreases for larger axis ratios. This is expected from the 
discussion at the end of Sect. 12.31 
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Fig. 8. Approximate equilibrium temperatures of 
spheroidal iron grains with a volume equal to a volume 
of a sphere with radius — 0.01 /im and = 0.05 /im 
as function of the ratio ajh. The solid and the dashed line 
correspond to calculations where the dependence of the di- 
electric function on grain size and temperature has been 
taken into account and has been ignored, respectively. 

At extreme values of the ratio ajh the equilibrium tem- 
peratures of spheroids with Vy ~ 0.05 /im are higher than 
for spheroids with — 0.01 /im. This can be explained by 
the higher collision rate due to surface scattering in the 
smaller grain leading to a higher electric dipole absorp- 
tion. 
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By including the size and temperature dependence of 
the dielectric function the temperatures of spheroids with 
Tv = 0.01 fim are up to '--^ 10 K lower than the tem- 
peratures where this dependence is ignored. If the depen- 
dences are not taken into account the temperature be- 
haviour would qualitatively be similar to the larger grains 
showing an increase of the temperature for intermediate 
values of the ratio a/b first and a decreasing temperature 
for more oblate and prolate shapes. 

3.1.3. Flickering small iron grains 

We derived the probability distribution of the tempera- 
tures of small grains using a numerical method described 
by Guhathakurta & Drainc (1989). To calculate the ther- 
mal energy of the iron spheres we made use of tabu- 
lated values of the heat capacity up to a temperature of 
298 K ('American Institute of Physics Handbook M972 ). 
At higher temperatures an analytical expression of the 
heat capacity is used (Chase EMJ- For this, we assumed 
iron to be in the a-5-phase. 

The temperature distributions of iron grains with ra- 
dius a = 10 A and a = 40 A heated by the ISRF (x = 1) 
are shown in Fig. El as dp(r)/dlog(T). For comparison 
we also included in this figure the theoretical tempera- 
ture variation of small spheres with a dielectric function 
of bulk iron at room temperature. If the size dependence 
of the collision rate is taken into account the probabilities 
of high temperatures are partly strongly reduced while the 
probabilities of low temperatures is clearly increased. Iron 
spheres with a = 10 A for example should have temper- 
atures lower than 10 K most of the time. In contrast a 
model where the size dependence of the dielectric func- 
tion is ignored predicts for the same grain size tempera- 
tures higher than 10 K with a most probable temperature 
of 30 K. 

The maximum temperature very small grains reach is 
roughly the temperature where the thermal energy is equal 
to the maximum photon energy of 13.6 eV. Because of the 
high heat capacity of iron grains the maximum tempera- 
ture is clearly lower in comparison with other grain species 
like silicate or graphite grains. By using the heat capac- 
ity given by Dwek |1986) graphite grains with a radius 
of 10 A attain a maximum temperature of above 400 K 
fFischera l2000|l . In contrast the maximum temperature of 
iron grains with the same size is less than 300 K. 

3.2. The SED of spherical iron grains in the ISM 

From the analysis of the extinction curve Mathis, Rumpl 
& Nordsieck H1977fl proposed a power law dn{a) cx a^'' da 
with k ~ 3.5 for the grain size distribution in the ISM. 
The minimum and maximum grain size as the power of 
the size distribution depends somewhat on the assumed 
grain composition. For the iron grains we considered a 
maximum grain size of Cmax = 0.1 /im. To study the effect 
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Fig. 9. Temperature distribution of spherical iron grains 
with radius 0.001 and 0.004 /im heated by the ISRF 
(MMP). Solid and dashed lines correspond to calculations 
in which we included and ignored the dependence of the 
dielectric function on temperature and grain size. 

of stochastic heated grains we assumed as minimum grain 
size different radii 10, 20 and 40 A. 

The mean intensity emitted at the surface of a flicker- 
ing grain with radius a is given by 

/,(a) = J dTp{a,T)B,{T)Q,^-,{a,iy,T), (8) 

where B^{T) is the Planck-function andp(a, T) the proba- 
bility that the grain with radius a has the temperature T . 
The spectral energy distribution (SED) vl^ of three flick- 
ering iron grains of different sizes is shown in Fig. 1101 By 
including the size and temperature dependence of the di- 
electric function the emission spectrum is shifted to longer 
wavelengths. This is expected as these grains cool more ef- 
ficiently and therefore, as shown in the previous section, 
have higher probabilities at low but lower probabilities at 
high temperatures. 

The total flux density of a distribution of spherical iron 
grains at distance D is given by: 

F^^^ j daAa-^"^ a?T:h{a), (9) 

where we normalized A to the total dust mass of 
iron grains with the density pd = 7.87 g/cm^: 

A^A'U/ j daa-^-'^^Tra^pd- (10) 

The SED vF^ of iron grains of 1 Mq at 1 kpc is shown 
in Fig. ^] For comparison wc give the SEDs for the cal- 
culation including and ignoring the size and temperature 
dependence of the dielectric function. Due to the large 
variation of the equilibrium temperatures the SED is rel- 
atively broad. In comparison to graphite or silicate grains 
the SED is less effected by flickering grains. A clear in- 
crease at the blue side of the spectrum requires grains 
with radii smaller than Omin ~ 20 A. The inclusion of the 
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Fig. 10. Emission spectra of individual small iron grains. 
The sizes are chosen to be 10, 20, and 40 A. The spectra 
shown as solid lines correspond to calculations in which 
the size and temperature dependence of the dielectric 
function has been taken into account. The spectra shown 
as dashed lines correspond to calculations in which this 
dependence has been ignored. 

size and temperature dependence of the dielectric func- 
tion reduces the emission at the blue side which is con- 
tributed by stochastically heated small grains. The peak of 
the emission on the other hand is shifted towards shorter 
wavelengths. At long wavelengths we obtained for the two 
different cases of the dielectric function almost the same 
emission spectrum. 



4. Discussion 

The study of the temperature behaviour of metallic grains 
in the ISM requires a realistic model of the ir optical prop- 
erties. The grains emit mainly in the IR, in which range 
the absorption is strongly affected by the free electrons. 
This leads to an absorption that depends on temperature 
and, in cases where the free path of the electrons is larger 
than the grain size, especially on the grain size. A simpli- 
fication can easily cause bigger uncertainties in the equi- 
librium temperatures of spherical and spheroidal grains or 
in the probability distribution of temperatures from flick- 
ering small grains. 

To derive the cooling behaviour of iron grains Tabak 
and Straitiff (IWT) and Tabak (HM?) only considered the 
temperature but no size dependence of the optical proper- 
ties. In addition it was assumed that the grain absorption 
is sufficiently described by the electric dipole approxima- 
tion. Both simplifications lead to absorption coefficients 
which are too low. Therefore, equilibrium temperatures in 
the range of 100 K for iron grains heated by the ISRF 
as presented by Tabak & Straitiff ( 1983J and also Tabak 
P987|) should be in general too high. 

The temperatures of iron spheroids with = 0.01 fim 
as given by Voshchinnicov et al. H1999|l show qualitatively 
the peculiar behaviour we derived only for larger spheriods 
where the absorption behaviour at long wavelengths of the 
corresponding sphere is dominated by the magnetic dipole 
absorption. The peculiar curve of spheroidal grains with 
0.01 fim should disappear if the size dependence of the 
dielectric function is taken into account. 




X [jum] 

Fig. 11. Spectral energy distribution (SEDs) of spherical 
iron grains with a mass of 1 Mq at a distance of Ikpc 
heated by the ISRF (MMP). The grains are assumed to 
have a power law dn(a) oc a"'^'^ da with a maximum grain 
size ainax = 0.1 /im and different minimum grain sizes Omin 
of 10, 20 and 40 A. SEDs in which calculation the temper- 
ature and size dependence of the dielectric function was 
included are shown as solid lines. Dashed lines correspond 
to SEDs in which this dependence has been ignored. 



4.1. Uncertainties in grain temperatures 

4.1.1. Uncertainties due to impurities and defects 

It is likely that the iron particles in astrophysical environ- 
ments show a high degree of impurities and lattice defects. 
And it is possible due to the surface tension that the lat- 
tice structure becomes more disordered towards smaller 
grains. Both impurities and defects lead to a higher scat- 
tering rate inside the solid in comparison to bulk iron. As 
shown this effect should be negligible for iron grains with 
radii smaller than 0.25 fim if C < 100. Larger grains on 
the other hand will be affected. 

We found by considering iron spheres heated by the 
ISRF that the grain temperature as a function of radius 
shows a minimum at ^ 0.3 /im and that the temperature 
increases towards even larger iron grains. This is caused 
by the fact, as is shown in Sect. 1X1 and Sect.0 that the 
absorption behaviour of big iron grains is dominated by 
the magnetic dipole absorption (Eq. IA.7|I which decreases 
towards larger grains if the scattering rate is dominated by 
scattering events at the grain surface and therefore given 
hy T — vp /a. In the case of less pure and disordered iron 
grains the collision rate will be dominated by scattering 
events at impurities and defects so that the absorption 
coefficient becomes essentially constant with grain size. It 
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is therefore possible that the increase in grain temperature 
will disappear if iron grains in astrophysical environments 
are considered. 

4.1.2. Uncertainties in the dielectric function 

A better description of the iron temperatures requires a 
better understanding of the optical properties in the IR. 
To include the temperature and size dependence of the 
dielectric function we assumed that the free electron part 
is given by a single Drude term. As a consequence we ob- 
tained high values in the IR which were attributed to the 
optical response of the bound electrons and were assumed 
to be independent of temperature and size. But measure- 
ments by Weaver et al. H1979|l of the optical properties of 
bulk iron at temperatures in the region of 140 K indicate 
that the contribution of bound electrons to the optical 
behaviour is strongly decreasing at energies lower than 
0.8 eV. 

In this case, as can be seen in Fig. |21 the IR-values of 
the dielectric function of very small grains could be lower 
and their absorption efficiency therefore higher than as- 
sumed here. If the part of the bounded electrons would be 
negligible so that the dielectric function is determined by 
the free electron part, the absorption or cooling efficiency 
of an iron grain with a — 0.001 /im could increase by al- 
most one order of magnitude. It is therefore possible that 
the probabilities at high temperatures of flickering small 
iron grains is even smaller than the probabilities presented 
here. 

5. Summary 

The optical properties of small spherical and spheroidal 
iron grains were derived using a self consistent model of 
the dielectric function including its dependence on size 
and temperature. We have shown that the limitation of 
the free path of the electrons due to the size of the grain 
has a non-negligible effect on the absorption behaviour of 
small spheres and spheroids. We applied the estimated ab- 
sorption coefficients for small spherical and spheroidal iron 
grains in the ISM which are heated by the ISRF and de- 
rived the emission spectrum of stochastically heated small 
iron spheres assuming as grain size distribution a power 
law dn (X a~^'^ da as proposed for grains in the ISM and 
chosing the radius of the largest grain as 0.1 /im. The main 
results for the temperature behaviour of small iron grains 
are: 

— The equilibrium temperatures of pure iron spheres 
heated by the ISRF are strongly size dependent. The 
highest equilibrium temperature of '--^ 50 K was found 
for iron spheres with radius ^ 0.01 /im and the low- 
est temperature of ~ 14 K for iron spheres with 
a ~ 0.3 /im. The temperature behaviour of larger iron 
grains is likely to be effected by impurities and de- 
fects. If these grains are relatively pure the tempera- 



ture should increase with size. Otherwise their temper- 
ature should be size independent. 

— In case of iron spheroids with a volume corresponding 
to a sphere with radius = 0.01 /im the temperature 
decreases for grains with oblate and prolate shape. 

By considering the size dependence of the dielectric 
function we found in comparison to calculations where 
this dependence is ignored: 

— The temperatures of small spheroidal grains {r^ ~ 
0.01 /im) should be significantly lower. 

— Absorption in the case of small grains (where electric 
dipole absorption is important) increases while absorp- 
tion in the case of bigger grains (where magnetic dipole 
absorption is important) decreases. 

— The absorption at resonance frequencies caused by 
shape effects of strongly elongated or flattened iron 
grains with = 0.01 //m is significantly reduced. The 
absorption at IR-wavelengths beyond the resonance 
wavelengths increases significantly. 

— Due to better cooling the probabilities of high temper- 
atures of small flickering grains should be lower and the 
probabilities of low temperatures significantly higher. 

— The contribution of flickering grains to the SED at 
short wavelengths is reduced. 
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Appendix A: Analytical approximation at long 
wavelengths 

Here we will discuss in some detail the absorption be- 
haviour of metallic grains using the electric and magnetic 
dipole approximation. To derive approximations it is as- 
sumed that the dielectric function is dominated by the 
free electron gas and therefore by . The cross section at 
wavelengths with 27ra/A <C 1 is given by (see e.g. Draine 
& Lee lTOMj) : 

Cab. = C:^s + C2s = ^Ini {a' + a™) , (A.l) 

where Im(a'^ -I- a™) is the imaginary part of the electric 
and the magnetic polarizability. In the case of spheres 
these are: 



- = --(^l-K-cot,--j w:th,^.e(-) (A.3) 
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If the contribution of the magnetic dipole absorption is 
negligible the absorption coefScient is known to be given 
by: 



Q 



abs 



Q 



abs 



CJjJt, 



(A.4) 



Therefore, the absorption is proportional to and in- 
dependent of grain size where the collision rate is given 
by r = vp/a. This is the case for very small iron grains 
at wavelengths larger ~ 60 fiia. At lower wavelengths the 
absorption is dominated by which causes a reduced ab- 
sorption with Qabs oc a/A^. 

With increasing grain size the magnetic dipole absorp- 
tion begins to dominate the absorption behaviour at long 
wavelengths. In the region |?/| ^ 1 and ^ F, so that 

~ iuip/ajT, the absorption coefficient due to induced 
eddy currents in the metal spheres can be simplified to 
the expression: 



i abs 



12 /Lua\^ ijJp 
90 \~) ' 



(A.5) 



Taking T = vp/o, the absorption is proportional to a'*/A^. 
Comparing the approximations of the electric (Eq. IA.4p 
and the magnetic dipole absorption (Eq. IA.5|) both terms 
are equal at a grain size of a « 0.015 /im. This results in 
the strong increase of the absorption at long wavelengths 
seen in Fig. 2|with grain size from 0.01 to 0.1 /xm. 

In the region <C F, where \y\ ^ 1, the absorption 
can be described by: 



Qabs ~ Q 



3\/2 



abs 



3 

£2 



3^2^^ 



VcjF 



(A.6) 



The absorption is therefore proportional to l/\f\ for a 
fixed collision rate F and proportional to l/\/aA if F = 
vp/a. For a pure cold iron sphere with a = 1 /im this 
behaviour occurs only at very long wavelength with A > 
1000 Aim (Fig.lBTj. 

The plateau with constant absorption coefficients at 
long wavelengths is caused by the magnetic dipole ab- 
sorption in the region |y| ^ 1 and ojp ^ oj ^ T where 
« —u)'^/u)^{l — iV /(jj). In this range the absorption co- 
efficients are approximately given by: 



i abs 



2c 



-V 



(A.7) 



Bigger iron grains in this parameter range absorb less ef- 
ficiently than smaller grains do if F = Wi^/a. As discussed 
in Sect. 14.1.11 this behaviour implies relatively pure and 
ordered iron grains with w 1. In case of less pure or 
disordered iron spheres the absorption efficiency may be 
independent of grain size. 

For completeness we also give the absorption be- 
haviour for the region tOp^ lo ^ T and \y\ <^ 1: 



Q 



abs 



12 w^a^F 
90 c3 



(A. 



The magnetic dipole absorption in this region is again 
independent of wavelength and is proportional to if the 
collision rate is given by F = vp/a. But this behaviour is 
not very prominent in case of iron spheres. Fig. 01 shows 
only a flattening of the absorption coefficients of a grain 
with a radius of 0.1 /im between 10 and 100 /im. 

Appendix B: Idealised dipole absorption of iron 

To visualise the effect of the limitation of the mean free 
path lengths of the free electrons on the cooling behaviour 
of iron spheres it might be useful to approximate the di- 
electric function by a single Drude term using the param- 
eters as given in Table El The derived absorption coef- 
ficients are shown in Fig. IB. II as function of wavelength 
and grain size. For comparison the absorption coefficients 
have been calculated for two different assumptions for the 
collision rate F. In the left hand figure no correction to 
the measured dielectric function (at room temperature) 
has been made. In the right hand figure it is assumed that 
the collision rate is dominated by surface scattering which 
may be appropriate to describe the situation of cold pure 
iron grains in the ISM which are heated by the ISRF. 

If no correction is applied, the magnetic dipole ab- 
sorption would be the dominant absorption at long wave- 
lengths for grains with radii larger than 3 nm. The limita- 
tion of the free path of the free electrons by the grain 
size shifts this value to grains with radii larger than 
^ 0.015 /xm. The absorption increases almost for all wave- 
lengths with grain size when no size dependence of the di- 
electric function has been considered (Fig. IB. II left hand 
side). When this dependence taken into account, the ab- 
sorption in the IR shows a maximum for grains with a 
radius in the range of ~ 0.3 /j,m (Fig. IB. II right hand 
side), which is the reason for the lowest temperatures of 
these grains in the ISM as derived in Sect. l3.lTn The de- 
crease of the absorption with grain size for grains larger 
than 0.3 /xm extends also into the region where the absorp- 
tion has to be described using the Mie-formalism (hatched 
region) . 

The approximations of the magnetic dipole absorption 
given in Sect. l2.2TTl are valid in the regions left of |y| w 1 
and right of \y\ « 3 (shown in Fig. IB. II as short dashed 
and dashed dotted lines) considering w 3> F and uo ^T. 

The absorption behaviour is shown for non-magnetic 
iron grains to demonstrate the different trends. As stated 
in the text, above A sa 1000 ^va the absorption of 
iron grains might be dominated by magnetisation effects 
(Fig. |SJ). At short wavelengths on the other hand the ab- 
sorption behaviour of the iron grains additionally may be 
effected by the response of the bound electrons. 
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Fig. B.l. Absorption coefficients QabslOjA) of metallic iron spheres derived approximating the dielectric function e 
through the free electron part e' . In the left hand figure the collision rate of the free electrons is taken to be the same 
as for bulk iron at 298 K. In the right hand figure it is assumed that the collision rate is dominated by scattering events 
at the grain surface. The absorption coefficients are shown as contours in logarithmic scale (thin solid lines) where the 
difference between the lines is taken to be A logj^Q((5abs(a, •^)) = 0.5. The labels correspond to log2Q(Qabs(fl, A)). The 
hatched region shows where the absorption behaviour deviates strongly from the dipole absorption (27ra/A > 0.5). The 
thick solid line visualises where the electric dipole absorption is equal to the magnetic dipole absorption. Where the 
frequency w of the photons is equal to the collision rate F of the free electrons in the grain is shown as long dashed line. 
The short dashed and dashed dotted line correspond to absorption coefficients where = 1 and \y\ =3, respectively. 
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